gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十四章/eg14_4.m

    %绘制系统的nyquist图
num=[20 20 10];
den=conv([1 1 0],[1 10]);
nyquist(num,den)

%进行局部放大
num=[20 20 10];
den=conv([1 1 0],[1 10]);
nyquist(num,den)                                        %绘制系统Nyquist曲线
axis([-2 2 -5 5])

%得到局部的nyquist曲线
num=[20 20 10];
den=conv([1 1 0],[1 10]);
w=0.1:0.1:100;                                             %指定频率范围
[re,im]=nyquist(num,den,w);                       %返回Nyquist曲线参数
plot(re,im)                                                    %使用Nyquist曲线参数绘制曲线
axis([-2 2 -5 5]);                                           %指定横纵坐标
grid;
title('系统(20s^2+20s+10)/[(s^2+s)(s+10)]Nyquist图(\omega>0)','fontsize',12);
xlabel('Real');ylabel('Imag')